Josephson 7r-state in a ferromagnetic insulator 
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We predict anomalous atomic-scale 0-tv transitions in a Josephson junction with a ferromagnetic- 
insulator (FI) barrier. The ground state of such junction alternates between 0- and 7r-states when 
thickness of FI is increasing by a single atomic layer. We find that the mechanism of the 0-7T transition 
can be attributed to thickness-dependent phase-shifts between the wave numbers of electrons and 
holes in FI. Based on these results, we show that stable 7r-state can be realized in junctions based 
on high-Tc superconductors with I^BaCuOs barrier. 
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PACS numbers: 74.50. +r, 72.26.-b, 85.75.-d, 03.67.Lx 



The developing field of superconducting spintronics 
comprises a plenty of fascinating phenomena that may 
complement nonsupcrconducting spintronics devices 
Mesoscopic hybrid structures consisting of superconduct- 
ing and magnetic materials have attracted considerable 
attention as devices with novel functionalities 0. One 
of most interesting effects is the formation of 7r-states 
in superconductor /ferromagnetic- metal/superconductor 
(S/FM/S) Josephson junctions jjjj]. Under appropriate 
conditions a fcrromagnet can become a 7r-phase shifter, 
providing the phase difference (f> = tt between two su- 
perconductors in the ground state in contrast to <t> = 
in ordinary Josephson junctions. Recently a quiet qubit 
based on S/FM/S 7r-junction [i| has been suggested as 
a promising device to realize quantum computation be- 
cause the spontaneously generated two-level system in 
this structure is robust against decoherence due to ex- 
ternal fluctuations. However, low energy quasiparticle 
excitations in a FM provide strong dissipation [f| . There- 
fore Josephson n junctions with a nonmctallic interlayers 
are highly desired for qubit applications @. Moreover, 
from the fundamental view point, the Josephson trans- 
port through a ferromagnetic insulator (FI) has been 
studied based on phenomenological models [7[ and not 
yet been explored explicitly. 

In this Letter, we study theoretically the 
Josephson effect in superconductor/fcrromagnctic- 
insulator/superconductor (S/FI/S) junctions using the 
tight-binding model. We show that the ground state in 
such structures alternates between the 0- and 7r-states 
when the thickness of a FI (Lp) is increasing by a single 
atomic layer. This remarkable effect originates from the 
characteristic band structure of a FI. Quasiparticles in 
the electron and hole branches acquire different phase 
shifts while propagating across a FI. We will show that 
the phase difference is exactly ttLf due to the band 
structure of a FI, thus providing the atomic-scale 0-7T 



transition. This mechanism is in striking contrast to 
the proximity induced Q-ir transition in conventional 
S/FM/S junctions. On the basis of the obtained results, 
we predict a stable 7r-state in a Josephson junction based 
on high-T c superconductors with a L^BaCuOs barrier, 
where electric current flows along the c axis of cuprates. 

Let us first consider an S/FI/S junction in the two- 
dimensional tight-binding model as shown in Fig. 1(a). 
The vector r = jx + my points to a lattice site, where 
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FIG. 1: (color online) (a) The Josephson junction with a 
ferromagnetic-insulator (FI) barrier on the tight-binding lat- 
tice. The magnetic moment in FI is chosen along the z axis 
in spin space. The band structure of a FI (b) and of standard 
band insulator (c) in the Bogoliubov-de Gennes picture. The 
dispersion for a hole with spin a is obtained as a mirror image 
of the dispersion for an electron with spin a with respect to 
Fermi energy. 
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x and y are unit vectors in the x and y directions, re- 
spectively. The thickness of a FI layer is Lp. In the 
y direction, we apply the hard wall boundary condi- 
tion for the number of lattice sites being W. Electronic 
states in a superconductor are described by the mean- 
field Hamiltonian, H = (1/2) 



g (c r h r r ' c r , 

(V2)Eres(4 A cl - T r A* cr) with 
c r = (crrf, c r ^), where c\ a (c r CT ) is the creation (an- 
nihilation) operator of an electron at r with spin cr(=^ 
or 4- )j c means the transpose of c, and <7o is 2 x 2 unit 
matrix. We introduce the hopping integral t among near- 
est neighbor sites and measure the length in the units of 
the lattice constant a. In superconductors, the Hamil- 



tonian leads h r 



-tSu 



+ (-fj, s + 4t)6 r ,r>}vo, 



the chemical potential /i s is measured from the band 
bottom and A = iAer 2 , where A is the amplitude of 
the pair potential in the s-wave symmetry and &j for 
j = 1 — 3 are Pauli matrices. We describe a FI by 
h r y = —t5u — r /i i&o ~ (.9/2 + 4t)8 rtr /(T3, where g cor- 
responds to a gap of a FI as shown in Fig. 1(b). 

The Hamiltonian is diagonalized by the Bogoliubov 
transformation. The Andrccv bound state consists of 
subgap states whose wave functions decay far from the 
junction interface. In what follows, we focus on the sub- 
space for spin-f electron and spin-J, hole [the dispersions 
shown by solid curves in Fig. 1(b)]. In superconductors, 
the wave function of a bound state is given by 
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where v = L (R) indicates a superconductor in the 
left (right) hand side, <p v is the phase of a super- 
conductor, = diag(e i ^/ 2 ,e" i ^/ 2 ), u{v) = [(1 + 
(-)VE 2 - A 2 /E)/2} 1 / 2 , and A, B, C and D are ampli- 
tudes of the wave function for an outgoing quasiparti- 
cle. The wave function in the y direction is xi( m ) = 
y/2/W sin[/m7r/(Fy + 1)], where / indicates a transport 
channel. The energy E is measured from the Fermi 
energy and k = cos _1 [2 — /i s /2i — cos{lir/(W + 1)} + 
i\J A 2 — E 2 / E] is the complex wave number. These wave 
functions decay as e~(J~ LF )/?° for j > Lp and eJ^° for 
j < with £0 being the coherence length. In a FI, the 
wave function is given by 
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where cosh/3 t = 1 + £7/2* + g/At + cos[lir/(W + 1)] - 
cos[W7r/(W+l)], cosh/?; = l+E/2t+g/At-cos[lir/(W+ 
1)] - cos[tt/(F7 + 1)], and JlJe^l and g R arc am- 
plitudes of wave function in a FI. The Andreev levels 



£n,i(4> — <t>L — 4>r) [n — 1, ••■ ,4] can be calculated 
from boundary conditions ^L(A,m) = ^ pi(X 1 m) and 
^r{L f + X,m) = ^>fi{Lf + A, m) for A = and 1. 
The Joscphson current is related to e n 1 via Ij(4>) = 

(2e/»)E»,i[^n,iW/^]/(en,i(^)). where f( £ ) is the 
Fermi-Dirac distribution function. The Josephson crit- 
ical current Iq is defined by Ic = Ij{' k /2). 

In Fig.[2)(a), we first show the Andreev levels e n .i = £ n 
for odd Lp (= 3 and 5) and even Lp (= 4 and 6) with 
W = I, n s = 2t, and A = O.Olt. The results show 
that the ground state for odd Lp is at (f> = tt, whereas 
that for even Lp is at <j> = 0. This atomic-scale 0-7T 
transition persists even if we increase Lp and W. In 
Fig. 2(b), we show the Josephson critical current as a 
function of Lp for W = 1. Temperature T is set to be 
O.OlTc -C T c , where T c is the transition temperature of 
a superconductor. The 7r(0)-state is always more sta- 
ble than the 0(7r)-state when the thickness of FI is an 
odd(even) integer. The reason is as follows. At low tem- 
peratures, only the Andreev levels below the Fermi en- 
ergy i.e., E\ and £2, contribute to Ic [see Fig. 2(a)]. 
In the odd (even) Lp cases, the 7T- (0-) state is stable 
because of d <j> £ 1 \ < p^/ 2 > (<)0, d^e 2 \^/2 < (>)0, and 
I d$E2 1 0=7r/2| > |^ei |^=W2 1- The atomic-scale 0-ir tran- 
sition is insensitive to W and material parameters such 
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FIG. 2: (color online) (a) The Andreev levels £i = Ei,i are 
plotted as functions of <f> for an one-dimensional S/FI/S junc- 
tion with W = 1. Left (right) panel shows the results for odd 
Lf (even L F ) with g = 0.5t. (b) Josephson critical current 
Ic at T = 0.01T C as a function of the thickness Lf for W = 1 
and g = 0.5£. The large red (small blue) circles indicate the 
7r(0) junction. In the inset, the 0-7r phase diagram on the 
g-Lp plane is shown for a two-dimensional S/FI/S junction 
with W = 10 at T = 0.01T C , where the red and blue regimes 
correspond to the 7r- and 0-states, respectively. 
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as (j, s , g, and A. As an example, inset of Fig. 2(b) shows 
the phase diagram on the g-Lp plane for W = 10. 

The mechanism of the 0-tt transition in a FI is very 
different from that in a FM. The key feature is expressed 
by the wave number of a quasiparticle in a FI as shown in 
Eqs. (UJ and (O, where q e and qh are the wave numbers 
for an electron spin-f and a hole spin-i, respectively. The 
real parts of q e and qi, reflect the wave number at the q 
points, where energy is closest to the Fermi energy, and 
differ by tt from each other . As shown in Fig. 1(b), the 
real part of q e is tt because the top of the electron band 
is located at q = tt. On the other hand, the real part of 
qh is because the top of the hole band is at q = 0. This 
is the origin of the difference between q e and qh which 
accounts the atomic-scale 0-tt transition. When we con- 
sider a usual band insulator as shown in Fig. 1(c), we 
always obtain q e = qh and their real parts equal tt be- 
cause both the top of the electron band and the bottom 
of the hole band are located at q = tt. As a consequence, 
0-state is always stable in usual band insulators. Thus 
we conclude that the characteristic band structure of a 
FI is the origin of atomic-scale 0-tt transitions. These fea- 
tures basically remain unchanged even when we consider 
Joscphson junctions in higher dimensions. In such junc- 
tions, however, the appearance of 0-tt transitions depends 
on relative directions between the current and the crys- 
talline axis. We will address this issue below. It should 
be emphasized that peculiar results presented above can- 
not be described by the standard quasiclassical Green's 
function method [8[ where band structure structure far 
from the Fermi energy is is ignored. 

Let us reconsider atomic-scale 0-tt transitions from a 
different view point of quasiparticle transmission coeffi- 
cient. In the high barrier limit {g 3> t), the Josephson 
critical current is perturbatively given by Ic ccT^T^ 
Here Ify,) is a transmission coefficient of a FI for an 
electron with spin-'f (-J,). By using the transfer-matrix 
method Q , T a for one-dimensional junctions can be ob- 
tained analytically T a w ai F (p a t/g) LF . Here p-^j.) = 
— (+)1 and «i F is a spin-independent complex constant. 
We immediately find T t /T 4 = (-1) Lf ■ Thus the relative 
phase of T a between spin-f and spin-J, is tt (0) for the 
odd (even) number of Lp. As a consequence, the sign of 
Ic oc (— 1) Lf becomes negative for odd Lp and positive 
for even Lp. In other words, a FI acts as a 7r-phase-shiftcr 
for the spin-f electron for odd Lp. 

The transfer-matrix method in real space also enables 
us to extend the calculations to another magnetic mate- 
rials. Up to now, we have considered uniform magnetic 
moment in FI, which can be schematically expressed by 
S/tit2 • ■ ■ tip/S or S/4ii 2 ■ • • Ilf/S. The arrows tj and 
\.j indicate the z-axis magnetization at j. We can ex- 
tend the above simple analysis to the arbitrary magneti- 
zation configuration, e.g., a random alignment described 
by S/WI2U ■■■ tL F -2tL F -iiL F /S. In such junctions, 
we find Ic ~ Yli=x,L P T* A T iA = IL=i, Lf (-1) = 



where Tj iCT is the transmission coefficient of an FI layer at 
i. Therefore we obtain a noticeable result, i.e., the sign of 
Ic is independent of magnetization configurations and is 
negative (positive) for odd (even) Lp. The appearance of 
the atomic-scale 0-tt transition has been also predicted in 
S/antiferromagnetic-interlayer/S junctions [10j. In their 
theory, however, the antiferromagnetic configuration is 
found to be essential for the atomic-scale transition. On 
the other hand, we conclude that the magnetization sym- 
metry is not necessary and that the 7r-phase difference 
between T-j- and Tj. is an essential feature for the atomic- 
scale transition. Therefore our analysis provides more 
general view for the physics of the atomic scale 0-tt tran- 
sition. 

Finally, we show the possibility of the atomic-scale 0-tt 
transition in a three-dimensional junction using realistic 
materials. Here we focus on La 2 BaCu05 (LBCO) [ll[ 
which is an important representative FI in spintronics. 
According to a first-principle band calculation [l2j| . the 
bottom of the minority spin band is at the T point whose 
wave number is (k a ,ki,,k c ) = (0,0,0), where kj for j = 
a, b, and c is the wave number along j axis (see Fig. 6 



in Rcf. 12). The mirror image of the minority spin band 
with respect to the Fermi energy corresponds to the hole 
band with minority spin in the Bogoliubov-de Gcnncs 
picture. Thus the top of the minority spin hole band is 
at the r point. On the other hand, top of the majority 
spin band is at the Z point with (k a ,kb,k c ) = (0,0, tt). 
Thus we can predict that the tt state would be possible if 
one fabricates a Josephson junction along c axis as shown 
in Fig. 4(a). Note that it is impossible to realize the tt- 
state if current flows in the a&-plane. This is because 
wave numbers in a6-plane at the bottom of the minority 
spin band and those at the top of the majority spin band 
are given by the same wave number (k a , k^) = (0, 0) [l2j ]. 

From the perspectives of the S/FI interface match- 
ing and the high-temperature device-operation, the us- 
age of high-T c cuprate superconductors (HTSC), e.g., 
YBa 2 Cu307_,5 and La2- x Sr x Cu04 (LSCO) is desirable. 
Recent development of the pulsed laser deposition tech- 
nique enables layer-by-layer epitaxial-growth of oxide su- 
perlattices [Hj]. In order to show the possibility of tt- 
coupling in such realistic HTSC junctions, we have cal- 
culated the c-axis Josephson critical current Ic based on 
a three-dimensional tight-binding model with L a and L^ 
being the numbers of lattice sites in a and b directions 
[Fig. 3 (a)] . In the calculation we have taken into account 
the ci-wave order-parameter symmetry in HTSC, i.e., 
A = Arf(cos k x a — cos k y a) /2. The tight-binding parame- 
ters t and g have been determined by fitting to the first- 
principle band structure calculations along the line from 
r to Z point [12| . Figure 3 (c) shows the thickness Lp de- 
pendence of I c at T = 0.01T C for a LSCO/LBCO/LSCO 
junction with g/t = 20, A^/t = 0.6, and L a = L\> = 100. 
As expected, the atomic-scale 0-tt transitions can be re- 
alized in such oxide-based c-axis stack junctions. 
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FIG. 3: (color online) (a) Schematic picture of a c-axis stack 
high-r c superconductor/FI/high-T c superconductor Joseph- 
son junction and (b) a d-wave ring to detect the 7r-junction 
behavior experimentally, (c) The Josephson critical current 
Ic as a function of the FI thickness Lf at T = 0.01T C for 
a c-axis stack LSCO/LBCO/LSCO junction with g/t = 20, 
Ad/t = 0.6, and L a — Lb = 100. The large red (small blue) 
circles indicate the tt(0) junction. 



The experimental detection of the 7r-junction is possi- 
ble by using a superconducting ring which contains two 
Josephson junctions as shown in Fig. |3Jb). When both 
junctions are in 0- (or n-) state at the same time, the 
critical current of the ring reaches its maximum at zero 
external magnetic flux. On the other hand, the critical 
current reaches its minimum at zero magnetic flux when 
the state is stable in one junction and ir is stable in the 



other [14 1 



From the view point of qubit applications, it is impor- 
tant to note that the harmful influence of midgap An- 
dreev resonant states 15|j 16| and nodalquasiparticles due 
to the d-wave symmetry on the macroscopic quantum dy- 
namics in c-axis HTSC junctions [17| is found to be weak, 
both theoretically [l8| and experimentally [l9|. There- 
fore we conclude that HTSC/LBCO/HTSC Tr-junctions 
would be promising candidates as basic-elements for 
quiet qubits. 

In summary, we have studied the Josephson effect in 
S/FI/S junctions based on the tight-binding model. We 
predict the formation of the atomic-scale 0-tt transitions 
in such junctions. This result is insensitive to the mate- 
rial parameters such as the gap g of the FI and the su- 
perconducting gap A, indicating that it is a robust and 
universal phenomenon. Our findings suggest the way of 
realizing ideal quiet qubits which possess both the quiet- 
ness and the weak quasiparticlc-dissipation nature. 
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